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The distribution of magnetic induction in Meissner state with finite London penetration depth 
is analyzed for platelet samples of rectangular cross-section in a perpendicular magnetic field. The 
exact 2D numerical solution of the London equation is extended analytically to the realistic 3D case. 
Data obtained on Nb cylinders and foils as well as single crystals of YBCO and BSCCO are in a good 
agreement with the model. The results are particularly relevant for magnetic susceptibility, rf and 
microwave resonator measurements of the magnetic penetration depth in high-T c superconductors. 

PACS numbers: 74.25.Ha, 74.25.Nf 
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The temperature and field dependencies of the mag- 
netic penetration depth yield basic information about 
the microscopic pairing state of a supercondiictorEI as 
well as vortex static and dynamic behaviorcm. Since 
most high-T c superconductors are highly anisotropic, a 
measurement in which the applied magnetic field lies 
at an arbitrary angle relative to the conducting planes 
yields a Meissner response arising from both in-plane and 
inter-plane supercurrents. The corresponding penetra- 
tion depths X a b and A c can differ widely in their magni- 
tude and temperature dependence and it is desirable to 
separate the two contributions to the measured penetra- 
tion depth. To study A b one must resort to a configura- 
tion in which the applied field is normal to the conducting 
planes so as to generate only in-plane supercurrents. Un- 
fortunately, the London equations in this geometry can- 
not be solved analytically, making it difficult to reliably 
relate the experimental response (typically a frequency 
shift or change in magnetic susceptibility) to changes in 
Xab. Exact analytical solutions are known only for spe- 
cial geometries: an infinite bar or cylinder in longitudinal 
field, a cylinder in perpendicular field, a sphere, or a thin 
film. These solutions are not practical since most high-T c 
superconducting crystals are thin plates with aspect ra- 
tios typically ranging from 1 to 30. Brandt developed a 
general numerical method to calculate magnetic suscep- 
tibility for plates and discso but this method is difficult 
to apply in practice and the solutions are limited to two 
dimensions. 

In this paper we describe the numerical solution of the 
London equations in two dimensions for long slabs in a 
perpendicular field. The results are then extended ana- 
lytically to three dimensions. We first compare our calcu- 
lations in the limit of A = with SQUID measurements 



on cylindrical Nb samples of differing aspect ratica. We 
then compare our calculations for finite A with data from 
Nb foils and platelets of both BSCCO and YBCO high-ii 
superconductors, obtained by using an rf LC resonatoru. 
Using numerical results and analytical approximations 
we derive a formula which can be used to interpret fre- 
quency shift data obtained from rf and microwave res- 
onator experiments as well as sensitive magnetic suscep- 
tibility measurements. 
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FIG. 1. Right half, gray scale image of the magnetic field 
in and around the sample of d/w = 1/5 and X/d — 0.5. Black 
color represents 5 = 0. Left half, contour lines of the vector 
potential. 



Consider an isotropic superconducting slab of width 
2w in the ^-direction, thickness 2d in the y-direction, and 
infinite in z- direction. A uniform magnetic field Ho is 
applied along the y - direction. In this 2D geometry the 
vector potential is A = {0,0, A}, so that the magnetic 
field has only two components H = {dA/dy, —dA/dx, 0} 
and the London equation takes the form: AA— \~ 2 A = 0. 
Outside the sample AA — —A-Kj/c — and dA/dn is 
continuous along the sample boundary. Here n is the 
direction normal to the sample surface. A numerical 
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solution of this equation was obtained using the finite- 
element method on a triangular adaptive mesh using 
Gauss-Newton iterations scheme. The boundary condi- 
tions were chosen to obtain constant magnetic field far 
from the sample, i.e., A (x, y) = —Hqx for y >> d and x 
>> w. 

Figure [j] presents the distribution of the the mag- 
netic field in and around the sample with w/d = 5 and 
X/d = 0.5. The black color on a gray scale image cor- 
responds to |B| = 0. The left half of the sample shows 
contour lines of the vector potential. Figure ^| shows pro- 
files of the y-component of the magnetic field at different 
distances y from the sample middle plane. 




FIG. 2. Profiles of the y-component of the magnetic field 
(parallel to the external field) for the sample shown in Fig. |l| 
Inset corresponding profiles of the vector potential. 

The inset shows the corresponding profiles of the vec- 
tor potential, normalized by its value A°(x = w) in the 
absence of a sample (a uniform-field curve = x is 
shown by the dotted line). Using the London relation 
47rA 2 j = — cA and the definition of the magnetic moment 
M = (2c) _1 JrX jd 3 r we calculate numerically the sus- 
ceptibility per unit volume (unit of surface cross-section 
in 2D case): 
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It is easy to check that for an infinite slab of width 2w in 
parallel field, where A = — XHo sinh (x/A)/cosh(w/A), 
Eq.(|l]) results in a known expression similar to Eq.(|j) 
below (with N = and R = w). In finite geometry 
there will be a contribution to the total susceptibility 
from the currents flowing on top and bottom surfaces. 
These currents are due to shielding of the in-plane com- 
ponent of the magnetic field, H x = dA/dy, appearing 
due to demagnetization. Figure ^ shows profiles of H x 
on the sample surface, at y — d, calculated for three dif- 
ferent samples, w/d — 8, 5, and 2.5. An analytical form 
for the surface magnetic field is known only for elliptical 



samples. We find, however, that it can be mapped onto 
the flat surface, so that the distribution of H x is given 
by: 



H x = 



(2) 



where r = x/w and a 2 = 1 + (2d/w) 2 . This equa- 
tion is similar to that obtained for an ideal Meissner 
screeningtm. Solid lines in FigJH are the fits to Eq.(||). 
The agreement between numerical and analytical results 
is apparent. 
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FIG. 3. Distribution of the in- plane, H x , component of the 
magnetic field on the sample surface, y = d. Symbols show 
result of numerical calculation and solid lines are the fits to 
Eq.@. 

Next we find a simple analytical approximation to the 
exact numerical results by calculating the ratio of the vol- 
ume penetrated by the magnetic field to the total sam- 
ple volume. This procedure automatically takes into ac- 
count demagnetization and non-uniform distribution of 
the magnetic field along sample top and bottom faces. 
The exact calculation requires knowledge of A(x, y) in- 
side the sample or H(x, y) in a screened volume outside, 
proportional to w 2 . The penetrated volume is: 



A |ff, 
Ho 



ds 



(3) 



where integration is conducted over the sample surface 
in a 31? case or sample cross-section perimeter in a 2D 
case. Using Eq.(||) for magnetic field on top and bot- 
tom surfaces and assuming H s = Hq/{1 — N) on sides we 
obtain: 
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Here N is an effective demagnetization factor and R is the 
effective dimension. Both depend on the dimensionality 
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of the problem. As mentioned earlier, Eq.(Q) is similar 
to the well-known solution for the infinite slab of width 
2w in parallel field. In that case R = w and the effective 
demagnetizing factor N = 0. In a 3D case (2w x 2w slab, 
infinite in the z— direction), R = w/2 and N = 0. The 
tanh(i?/A) term in Eq. (^)was inserted to insure a cor- 
rect limit at A — * oo. This correction becomes relevant 
at X/R > 0.4, which is realized only at about T/T c > 0.9 
for typical high-T c samples. 




FIG. 4. Calculated — 4ty\ (X) for a slab of w/d = 5. Solid 
line is a fit to Eq.(^J) with the effective dimension R/w — 0.36. 
Dotted line is calculated using R/w — 0.39 from Eq.(^) and 
a dashed line is a plot with R/w = 1. Inset: — <inx (A — > 0) 
calculated for samples of different aspect ratio. Solid line is 
1 + w/d. 

For the actual geometry studied here, both R and N 
depend upon the aspect ratio w/d. Unlike the case of an 
elliptical cross-section, the magnetic field is not constant 
within the sample so there is no true demagnetizing fac- 
tor for a slab. However, N can still be defined in the limit 
of A -> 0, through the relation, 4irM/V s = -H/ (1 - N). 
We find numerically that in a 2D case, for not too large 
aspect ratio w/d, 1/(1 — N) w 1 + w/d. Calculating the 
expelled volume as described above, the effective dimen- 
sion R is given by: 

w 

Rid = — . , _ n (5) 

1 + arcsm (a 1 ) 

In the thin limit, d -C w (a — > I), we obtain i?2D ~ 
0.39w. 

The natural extension of this approach for the 3D disk 
of radius w and thickness 2c? leads to 1/(1— N) » l+w/2d 
and 

w 

R3D = 2(1 +[! + (£)'] arctan(^)-^) (6) 

In a thin limit, i?3u ~ 0.2w. Eq.(6) was derived for a 
disk but the more experimentally relevant geometry is a 
rectangular slab. There is no analytical solution for the 
slab. However, a 2 = 1 + (2d/w) 2 is relatively insensitive 
to w in the thin limit and so we approximate w for a slab 



by the geometric mean of its two lateral dimensions. The 
validity of this approach will be determined shortly. 

To verify Eqs.(@) and (|) we calculated \ W numeri- 
cally. The result is shown it in Fig. || by symbols. The 
solid line is a fit to Eq.(|) with N = 0.86 and R/w = 0.36. 
The effective dimension calculated using Eq. (||) gives 
R/w = 0.39 and the corresponding susceptibility curve 
is shown as a dotted line. The calculated effective de- 
magnetization factor is N — 0.84. It is seen that our 
approximations are reasonably good, ft should be borne 
in mind that these are all 2D results - the sample extends 
to infinity in the z-direction. Demagnetizing effects are 
significantly larger in two dimensions than in three owing 
to the much slower decay of fields as we move away from 
the sample (compare 3D sphere, N = 1/3, and cylinder 
in perpendicular field, N = 1/2). Therefore we expect 
our approximations to work better in three dimensions. 
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FIG. 5. Linear magnetic susceptibility of Nb cylinders of 
different aspect ratio measured at T = 8 K. Solid line is a 
plot of 1 +w/2d 

In a 3D case the validity of our results can be verified 
experimentally by independently measuring the demag- 
netization factor as a function of the aspect ratio and the 
magnetic susceptibility for a finite London penetration 
depth A. To achieve the first goal, we measured niobium 
cylinders of radius w and length 2d using a Quantum 
Design MPMS-5 SQUID magnetometer. Sample dimen- 
sions were typically of the order of millimeters, which 
allows us to disregard London penetration depth of Nb 
(about 500 A). The initial susceptibility obtained from 
the magnetization loops at T — 8 K is shown in Fig. ||. 
The solid line is a plot of l + w/2d (not the fit) and for an 
aspect ratio up to w/d = 10 the agreement is excellent. 

To test our result for R (Eq.||) in actual samples we 
need the magnetic penetration depth. It is common 
to measure changes in the penetration depth by using 
the frequency shift of a microwave cavity or an LC res- 
onator. In these techniques, the relative frequency shift 
(/ ~ /o) //o due to a superconducting sample is propor- 
tional to H~ 2 J M ac • HdV, which in turn is proportional 
to the sample linear magnetic susceptibility (M ac is the 
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ac component of the total magnetic moment, H is the 
external magnetic field and fo is the resonance frequency 
in the absence of a sample). Using Eq.([|) and Eq.(^|) we 
obtain for A << R: 



A/ _ V. 
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(7) 



where V s is the sample volume, Vq is the effective coil vol- 
ume. The apparatus and sample - dependent constant 
A/ = Vs/o/(2Vo (1 — N)) is measured directly remov- 
ing the sample from the coil. Thus, the change in A with 
respect to its value at low temperature is 
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FIG. 6. Effective dimension 7? determined experimentally 
for different samples (symbols). The upper solid line is an 
"infinite slab" model (R = w/2) and the lower solid line is an 
analytic approximation R ~ w/5. 

We used an rf tunnel-diode resonator! to measure Sf 
in Nd foils, YBCO and BSCCO single crystals. Com- 
bining Sf with an independent measurement of AA(T) 
and a measured value for A/o, we then arrived at an ex- 
perimental determination of the effective dimension R. 
For the Nb and YBCO samples, AA(T) was obtained 
using the demagnetization-free orientation (rf magnetic 
field along the sample ab— plane) where R — w and 
1/(1 — JV) = 1. In BSCCO, the large anisotropy pro- 
hibits using this method and we used reported values 
of dX/dT ~ 10 A/KBq. Figure Fig. || summarizes our 
experimental results. The upper line represents the "in- 
finite slab" model, where R = w/2, whereas the lower 
solid line is R = 0.2w obtained in a thin limit of Eq. 
(||). Symbols show the experimental data obtained on 
different samples, indicated on plot. In three samples: 
YBCOl (w/d = 57), Nbl (w/d = 29) and Nb2 (w/d = 
15), R agrees with Eq.(||) to better than 5 %. The stan- 
dard result, R = w/2, is too large by a factor of 2.5. 
Both YBC02 and BSCCO give R roughly 20 % smaller 



than predicted. For the BSCCO data, it is possible that 
a sample tilt combined with the very large anisotropy of 
A produces an additional contribution from A c . If the c- 
axis is tilted by an angle, 9 away from the field direction, 
the frequency shift is given by 
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The importance of the tilt depends upon the relative 
changes in X a b and A c with temperature. From Eq.(j^) 
we obtain for the relative contribution to the frequency 
shift: 
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where we used the previous estimates of N and R. For 
BSCCO we take, d\ c /dT ~ 170 A/K and dX ab /dT ~ 
10 A/KM, Eq. © reduces to « 1 + tan 2 (0). Wc 



then find that for sample tilt to produce an additional 
20 % frequency shift a misalignment of 8 20° would 
be required. Our estimated misalignment was a factor 
of 10 smaller than this so the discrepancy between mea- 
sured and predicted R was not due to tilt. Both the 
BSCCO and the YBC02 sample were more rectangular 
than square and our use of the geometric mean for w 
could be the source of the error. 

In conclusion, we solved numerically the London equa- 
tions for samples of rectangular cross-section in perpen- 
dicular magnetic field. We obtained approximate formu- 
lae to estimate finite-A magnetic susceptibility of platelet 
samples (typical shape of high-T c superconducting crys- 
tals). 
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